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расчётным способом. Метод основан на гипотезе, теоретическое обоснование кото- 
рой осуществлено с позиций теории информации. Рассмотрен метод формирования 
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Введение. Под математическим моделированием понимается такой способ 
моделирования, при котором математическая модель того или иного объ- 
екта (или системы объектов) представляется в виде некоторых аналитиче- 
ских зависимостей. Такой подход к моделированию позволяет получать 
результаты в наглядном виде, применять модели одних объектов (напри- 
мер, в сети городского пассажирского транспорта) при создании модели 
других объектов (органов планирования и управления) [1]. 

Математическое моделирование позволяет использовать богатый 
набор оптимизационных методов, доказать корректность применения этих 
методов. Кроме того, намного упрощается анализ получаемых результатов. 

Основным недостатком математического моделирования является 
ограничение на число факторов, которые можно учесть в модели. При 
большом числе факторов модель усложняется настолько, что теряется вся- 
кая наглядность получаемых результатов. 

Постановка задачи. При описании модели мы будем рассматривать зада- 
чи, решаемые на стадии текущего планирования и управления транспорт- 
ными процессами, а именно: 

1) выбор оптимальной системы маршрутов; 

2) оптимальное распределение подвижного состава по маршрутам; 

3) составление оптимального расписания и др. 

Основные методы. В данной статье рассматриваются методы формиро- 
вания информационного обеспечения перечисленных задач, а также пред- 
лагается метод решения задачи оптимального распределения подвижного 
состава по маршрутам. 

Для формирования информационного обеспечения необходимо рас- 

смотреть теорию расчетов пассажиропотоков, в том числе методы оценки 
существующих матриц корреспонденции [3]. 
Модель восстановления матрицы корреспонденций. Наиболее полно 
методы формирования и оценки матриц корреспонденций рассмотрены в 
отчете ПАЗА [5]. Методы оценки существующих матриц корреспонденций 
предлагается классифицировать следующим образом: 
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- методы, основанные на моделировании распределения транс- 
портных потоков (Тга#Ис тоде!тд Базе арргоасйе5); 

- статистические методы (З4ай$Иса! [\еегепсе арргоаспе5); 

- градиентный метод (СгаЧепе Базе зоиоп {есйпюче5). 

Автор этого отчета Т. АбгаПатзоп считает, что в наиболее общем 
виде модель обновления матрицы корреспонденций формулируется как 


ттЕ(8,>)=у,Е(5,8)+ УЕ, (у), (1) 
у, >0,у= а551еп(е), 
где © - новая матрица корреспонденций; © -— матрица корреспонденций, 


полученная по результатам предшествующих исследований; 
иу - вектора значений наблюдаемых потоков — новый и прежний; 
Ет, Е›- некоторые меры расстояния, например, евклидова норма, или 
минимум информации (максимум энтропии); а551еп(2) — распределе- 
ние потоков у по сети, разделяющее матрицу корреспонденций г, на 
потоки, следующие разными маршрутами. 

В случаях, когда отсутствует “старая” матрица, используют методы, 
основанные на моделировании распределения потоков. В этих методах 
наиболее сложным является установление маршрутов, по которым следуют 
корреспонденции д, из зоны /в зону ]. Такая информация представляется в 


виде матрицы Р, элементы ру которой являются долей корреспонденций 


между зоной /и зоной 7, использующих дугу а 


[0 
Оз р; <1. (2) 
Соответственно интенсивность движения у, на данной дуге А 
__ а 
=» рус,’ @ ©А. (3) 
И 


В задачах восстановления матриц корреспонденций выделяют два 
типа распределения потоков 
1. Пропорциональное распределение. В этом случае потоки у, и 


доли ру считаются независимыми друг от друга. Значения Ру могут за- 


даваться до начала моделирования распределения. В том числе допускает- 
ся применение принципа “все или ничего” (а|-ог-по та), при этом пред- 
полагается, что корреспонденции используют пути минимальной стоимо- 
СТИ. 

2. Равновесное распределение. Моделирование выполняется в со- 


ответствии с принципами Водропа [4], т.е. значения Ру зависят от вели- 


чин потоков на всех дугах р’; = ри(у) и не могут определяться заранее. 


Отсутствие предварительной информации о маршрутах движения 
требует представить исходные данные для оценки матрицы корреспонден- 
ций так, чтобы каждая пара взаимно корреспондирующих пунктов была 
связана одним маршрутом. Это можно осуществить разделением матрицы 
корреспонденций участка дорожной сети на группу матриц корреспонден- 
ций, каждая из которых будет представлять “отдельный маршрут”. Анало- 
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гом является модель оценки межостановочных корреспонденций отдельно- 
го маршрута. 

В модели изолированного маршрута рассматривается одно из двух 
направлений движения и используются данные о количестве входящих и 
выходящих пассажиров на остановочных пунктах (ОП); межостановочная 
матрица корреспонденций Х является треугольной с нулевой диагональю, 


0 дв в зе 
О х И" 
23 о 2,п > 
(4) 
0 Хип? 
0, 


где х, - количество пассажиров, проехавших между / и / остановочными 
пунктами, х, >0. 


Сумма элементов #й строки матрицы соответствует количеству пас- 
сажиров, вошедших на /ом остановочном пункте, а сумма элементов го 
столбца матрицы соответствует количеству пассажиров, вышедших на ом 
остановочном пункте 


а=Ух, т 1 =1|,....П, (5) 
1 


| 
при этом а;и Б; удовлетворяют условию 


У а -Уь. (6) 


Для матрицы Х, максимизирующей взвешенную энтропию и отве- 
чающей условию 





А. 
Н(х)= Ух, - > тах, (7) 
й Х; 
было получено решение в следующем виде: 
Х!2 = Б?, 
ры 
а, — 2х Ь, 
х, = НЫ —. 0: (8) 
Е 
Л 
() ‚ еси С. =0; 
а_Б. 
_ "17 если С. >00: " 
() ‚ если С, =0; 
= Я | 
где С — — Чт — р та’. 
КЕ и=К+1 


Определенные положения рассмотренной модели оценки межоста- 
новочной матрицы корреспонденций изолированного маршрута могут быть 


514 


Вестник ДГТУ, 2009. Т9. №3(42) 


использованы в модели восстановления матрицы корреспонденций транс- 
портных потоков. 

Проведение эксперимента. Апробация модели проходила на предпри- 
ятии МУП ПО «Электроавтотранс» г. Усолье-Сибирское. В городе существу- 
ет три маршрута: №1 - 11.9 км, №2 - 18,5 км, №3 - 12,8 км. Маршрут №2 
в межпиковые часы движется по укороченному маршруту. В будние дни на 
линию выходят 32 вагона, в выходные - 28. Интервал движения составляет 
5 минут в будние дни. 

В силу относительной простоты и малой трудоёмкости способ на- 
блюдений сегодня представляется наиболее практичным. Он даёт возмож- 
ность многократных повторов и быстрой обработки результатов. Обследо- 
вания этого типа не требуют вхождения в контакт с пассажирами. Практи- 
ковались следующие наблюдения: 

- подсчёт посадок и высадок пассажиров на остановочных пунктах 
учётчиками, находящимися в подвижном составе; 

- глазомерная оценка водителями, кондукторами или специальными 
учетчиками наполнения подвижного состава на перегонах (в сечениях). 

В условиях города Усолье-Сибирского с учетом указанных оценок и 
реальных возможностей было решено использовать наблюдение. Наблюде- 
ние с целью глазомерной оценки наполнения подвижного состава можно 
производить внутри подвижного состава или извне. Целесообразно исполь- 
зовать сочетание этих способов. Мы использовали наблюдение внутри под- 
вижного состава. В результате была получена следующая матрица пасса- 
жирских корреспонденций (табл. 1). 

Таблица 1 
Элементы межостановочной матрицы пассажирских перевозок 


Номер оста- Номер остановки выхода 
новки входа 





В табл.1: а; — число пассажиров, вошедших в трамвай на ЁиИ оста- 
новке; 6; — число пассажиров, вышедших из трамвая на }й остановке; 
х; - элементы межостановочной матрицы пассажирских корреспонденций, 
показывающие число пассажиров, проехавших от #й до Ги остановок. 

Тогда сумма всех элементов в #й строке соответствует числу пас- 
сажиров, вошедших в трамвай на /й остановке (а), сумма всех элементов в 
ом столбце соответствует числу пассажиров, вышедших из трамвая на 
том остановочном пункте (Ъ)}. 

Визуальный метод позволяет вначале получить элементы матрицы 
корреспонденций пассажиров, а затем, суммируя их по строчкам (столб- 
цам), получаем величины вошедших (вышедших) пассажиров на каждом 
остановочном пункте всего маршрута (табл.?2). 
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Таблица 2 
Сводная таблица обследования пассажиропотока на маршруте №1 
(прямое направление) 


Название остановки выхода 


Название ос- 
тановки входа 


оборудования 


Трампарк 
(кольцо) 
Загородная 
Кристалл 
Юбилейный 
Гостиница 


ЕЕ Бан 


Не Магистральная 


Е 

ЕЕ 

М О О О 
Ленина ВЕ ОИ ПОВ ЕВЕ 
Институт __ МИНИ 

Магистральная _5 [ПГ Пин 


База оборудо- 3 
вания 


9-я столовая 1 | 
Кю Ор м 
[Общее кол-во | 115 [| ПИ 


Схожие данные получаются и при обследовании пассажиропотока 
на обратном направлении. 

По данным о количестве вошедших и вышедших пассажиров на ка- 
ждой остановке всего маршрута был предложен метод определения пасса- 
жирских корреспонденций расчётным способом. Метод основан на гипоте- 
зе, теоретическое обоснование которой осуществлено с позиций теории 
информации. 

Пусть Р,- вероятность того, что пассажиры, вошедшие в трамвай 
на ий остановке, выйдут на остановке } т.е. 


| ч 
= 5 
= |9. |5 
|= 
-- о ь 4 
36| 
_22 | 9 
18| 
И И и 
Е а 
Ще 


= 

-П 

т 

25 

Ф 

= 

и ®) 

х 
ЕЕ 
аа 
ВН 
ть 


[6% 
Г. 
> 
т 
ыч 
ыы 
ИИ 
"ТЕТ 
ВИН 
ВЯ 





Ш „а 

Р.=1,  Р, >0,1=1п1-1 (10) 

И у 

1=-4 
где п - число остановок на маршруте. 
Ш 
Предполагая, что х, — попарно независимы и 3 - =а,, можно 
= 
определить энтропию распределения вероятностей Р; относительно /- 
остановки 


Н‚ =—а, УР ‚108 Р, (11) 


1=НЫ 
Суммарная же энтропия равна 


Н.Н, = Ур 08,Р.; (12) 


=+/ 
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Р, — неизвестны, тогда максимизация суммарной энтропии определяет не- 
которое множество 


п-1 
Опар ; УР ве Ти- ЛЕЙ. 50) 
Я У=+1 
которое задаёт такую структуру искомой матрицы корреспонденций, когда 
известна лишь информация о количестве вошедших и вышедших пассажи- 
ров на каждой остановке всего маршрута. 
Так как х,= а,Р;, то 


— В» Ух р в, 7 > тах (14) 


1=| ]=7+1 
при ограничениях (5), (6). 
Предложенная методика основана на решении системы линейных 
уравнений вида 


= 
Ч; 5; 
п-1 
1=1 
УЗ: Ь 
К=-+1 
где [@ — некоторые неотрицательные числа. При этом элементы матрицы 
пассажирских корреспонденций определяются следующим образом: 


хе ри, ЕРЬ, (16) 


й 
о, 


С учетом этих о НобЗВЫЯ построены следующие рекуррентные 
соотношения: 


=5,, У=Т+Ь,П, (15) 








а. Хы: 0, 
1 
х, = НЕ к ЕСИ: 60. (17) 
. 
0 ‚ если С, =0; 
а.6, 
ОВ 8 ба. а 
0 ‚ еси С. =0; 
1—1 
где С, = ›3 Я т а 
К<1-—1 т=к+1 


Если предположить, что выход на /-Й остановке любого пассажира 
из перевозочной единицы из всех подъехавших ((;) есть независимое со- 


1 
бытие, вероятность которого равна — , то для определения пассажирских 
р 
корреспонденций расчётным способом можно предложить следующие 
формулы: 
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а) условное математическое ожидание случайных величин 4, при 
1—1 
выполнении условия а, =а, — 5 - усреднение по вероятности: 


т=1-+1 у вы 
=м(я, =, (19) 
К Со 
р Ь, —КА, если а, > Ь,; ое -4] если а, >в: 
где Я, = р 
а, — К, если а, <Б., К= 
й я а, — пи а,,(О, -Ь, а 
С... - Ул 


у 
Б) среднее арифметическое значение числа пассажиров, имеющих 
возможность проехать с некоторой вероятностью между двумя любыми 
остановками 
х,= т (20) 
Е =Ь, - шт ,(О,-а, | К =Ь,, если а, >Ь; 
К =а, — пил [а, (О, -Ь,)| К, =а,, если а, <Ь;; 
С) те значения К, вероятность которых наибольшая: 
Е 6, 
С. ем 
Ь. 
Со, 
где К определяется так же, как в пункте а). 
Полученные формулы (14) - (21) позволяют определить величины 
х,. По этим данным проводились экспериментальные расчёты, в результате 
которых выяснилось, что отклонение численных значений колеблется в 
достаточных пределах, удовлетворительных для большинства практиче- 
ских расчётов. 
Программное обеспечение. По результатам эксперимента создана про- 
грамма, которая позволяет рассчитать потребное количество трамваев, их 
марку и вместимость для заданного пассажиропотока как на маршруте в 
целом (максимальное количество транспорта), так и по часам суток (по- 
требное количество транспорта на каждый час). Исходными данными для 
расчета являются: 
а) максимальный объем перевозки пассажиров в час «пик», пас/ч; 
6) суточный объем перевозки пассажиров, пас; 
в) количество часов работы трамвая на данном маршруте за су- 
тки, ч; 
г) время оборота, ч. 
Основным критерием при выборе типа транспортного средства для 
перевозки пассажиров на заданном маршруте является целесообразный 
интервал движения, который определяется по результатам анализа пасса- 


где 


х, = тах (21) 
[4 
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жиропотока. Кроме интервала движения, учитывается номинальная вме- 
стимость подвижного состава, которая соответствует числу мест для сидя- 
щих пассажиров и числу мест для стоящих из расчета 5 человек на 1м* 
свободной площади салона. При выборе типа трамвая (по рациональному 
использованию пассажировместимости) задается три значения интервала, 
соответствующие рангу маршрута (табл.3), прибавляется одна минута и 
определяется вместимость. 

Таблица 3 

Ранжирование маршрута 


ПО Пани; Е Свыше 3000 
935 1500...2500 


ош бб 600...1000 
11...20 300...500 
21...80 До 300 


Затем отнимается одна минута и расчет повторяется. Расчетная 
пассажировместимость сравнивается с общей вместимостью (см.табл.3), и 
выбирается транспорт, вместимость которого будет использоваться более 
рационально. А, следовательно, будет наблюдаться снижение себестоимо- 
сти перевозок. 

В результате обработки данных о маршруте на экран выводятся 
(рис.1): марка транспорта; вместимость, пас; максимальное количество 
транспорта для часа «пик». 





`“ Определе ние потре бности в подвижном составе 


Файл Вычисления 
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Рис.1. Определение потребности в подвижном составе 
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При работе с программным обеспечением в справочные таблицы 
вводятся все марки транспорта, в частности значения Б, которые сущест- 
вуют на автотранспортном предприятии (рис.2), а также таблицы коэффи- 
циентов а; — суммарное количество пассажиров, перевозимое всеми транс- 
портными средствами на каждом маршруте для каждого временного интер- 
вала. Кроме параметров а;, для каждого временного интервала указывают- 
ся значения ар, $; (рис.3). Значения $; - это оценки стандартов величин а; 
(точности их определения), они требуют для оценки качества решения - 
точностей определения параметров модели коэффициентов корреляции 
между ними. Коэффициенты характеризуют взаимовлияние параметров - 
возможность увеличения одного за счет уменьшения другого. 


Инте ;звалы вренени 









Код 


_ Промежутки времени, ч 
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Рис.2. Лист НЫ таблицы» 





| ‚ Начиная со строки л. Статистическое | 
|! Коэффициенты а; для временного интервала, максимальные количества 0! р ы те 
моделирование для оценки качества | 
‚| перевозимых пассажиров на маршрутах и соответствующие т и 2: нс 
а Е : 
и 





Рис.3. Исходные параметры модели ах», ду, $; (1 =1,2,3, ]=1,2,3) 
При работе с программой надо перейти на рабочий лист (рис.4), 


выбрать временной интервал и выполнить вычисления, нажав кнопку «Вы- 
полнить вычисления». 
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Здание а Не а аа И иене феями фуцжьнненияке 


Расчет оптимального числа транспорта типа } (=1.. п) на маршруте 1 (1=1..ла) для заданного вр 


Выбор интервала 
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т те лист ит инт. 1 у Инт. _2 Ки инт. 3 а Е мнт. — ву 





Рис.4. Рабочий лист 


Программа сама подберет, учитывая заданные параметры, необхо- 
димое количество транспортных единиц на каждый маршрут. Сравнивая 
два интервала: 1-й и 2-й (см.рис.2), получаем сравнительные таблицы рас- 
четов (рис..5). 










ян + ен Здон-и не --—-- 


Выбор интервала 


‚ Выбор интервала 
о Выполнить вычислени 


= тит 2 


Значения ХВ - см. постановку задачи 





Значения Х\ - см. постановку задачи 


анны 


Рис.5. Сравнительные таблицы расчетов на разных временных интервалах 


Из рис.5 видим, что в часы пик количество транспортных средств 
значительно увеличивается, а в первый временной интервал можно не вы- 
пускать на линию трамвай марки 71-619 КТ; это хорошо видно и на диа- 
грамме, которая построена по результатам исследований (рис.6). 
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Количество транспортных средств 
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Рис.6. Результаты исследований 


Выводы. С помощью разработанной программы определения потребности 
подвижного состава можно судить об оптимальности парка подвижного 
состава города и возможности прогнозирования требуемого парка подвиж- 
ного состава города с точки зрения технико-эксплуатационных показате- 
лей, о выборе оптимальной системы маршрутов, оптимального распреде- 
ления подвижного состава по маршрутам и о составлении оптимального 
расписания. 
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